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1. Introduction 

An explanation of the origin of intermittency in the fully developed turbu- 
lence remains one of the main open problems of theoretical hydrodynamics. 
Quantitatively, intermittency is measured by deviations from Kolmogorov's 
scaling of the velocity correlators. Below, we shall describe a recent progress 
in the understanding of intermittency in a simple system describing the pas- 
sive advection of a scalar quantity (temperature, density of a pollutant) by 
a random velocity field. The system, maintained in a stationary state by a 
large scale source, exhibits a down-scale energy cascade. If the molecular 
diffusion is small, an intermittent inertial range sets in [1]. We shall explain 
this phenomenon by identifying the origin of the anomalous scaling of scalar 
correlators in a simple model of passive advection due to Kraichnan [2]. 

2. Kraichnan model 

The advection of a scalar quantity 9(t, r) by a velocity field v(i, r) is de- 
scribed by the linear differential equation 



where f(t,r) represents the external source of the scalar. If the source 
vanishes then the scalar is carried along by the flow: 
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where r(t) is the (Lagrangian) trajectory of the fluid particle located at 
time to at ro: 

dr 

-=v(t,r), r(t )=r . (3) 

The non-zero source / keeps creating the scalar along the Lagrangian tra- 
jectory and Eq. (2) is modified to 

9{t, r(t)) = 9(t , r ) + ids f(s, r(s)) . (4) 

In the presence of diffusion of the scalar the above equations require small 
changes: 

1. Eq. (1) picks up a term kV 2 6 on the right hand side with k standing 
for the diffusion constant, 

2. Brownian motions (3{t) should be superimposed on the Lagrangian tra- 
jectories by adding the term to the velocity in Eq. (3), 

3. the right hand sides of Eqs. (2) and (4) should be averaged over (3. 

The Kraichnan model of the passive advection of the scalar [2] assumes 
velocities decorrelated at different times, and, for a fixed time, distributed 
as a Gaussian (i.e. non-intermittent) field with mean zero and a non-smooth 
typical behavior in space 

|v(t,r) -v(i,r')| - |r-r'|«/ 2 (5) 

where £ is a fixed parameter between and 2. This may be achieved by 
imposing the velocity covariance 

(v Q (i,r) v /3 (t',r / )) = S(t-t') V a P{r-r') (6) 

with 



p°/3( ) _ p^( r ) = D r f (d - 1 + £) <5 Q/3 - £ 



= D a ?(r) (7) 



where d denotes the space dimension. The typical behavior (5) follows since 
D a P(r) (x r^. The incompressibility V • v = of the velocity field is assured 
by the relation 9 r «D a "(r) = 0. 



We shall be interested in the effect caused by a steady injection of the 
scalar on distances longer than some large scale L. One may conveniently 
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model such a source by a random Gaussian field f(t,r) with mean zero and 
covariance 

</(i,r) /(*',!■')) = 8{t-t')C{^f) (8) 

where C(^) is approximately constant for r <^ L and decays fast for r>L. 
We shall assume / independent of the velocities v. 

3. Steady state 

The evolution equations of the hydrodynamical type imply identities for the 
correlators of the evolving quantities known under the name of Hopf equa- 
tions. These equation usually couple the correlators with different number 
of points. The case of the passive advection is no exception since its Hopf 
equations couple the scalar equal-time correlators 

N 

F N (t;r) = (l[8(t,r n )) (9) 

n=\ 

to the mixed correlators involving both the scalar and the velocity field. 
However, in the Kraichnan model with temporarily decorrelated velocities, 
the mixed correlators may be expressed by the known 2-point function of 
the velocities and the correlators of the scalar alone. The resulting Hopf 
equations may be easily seen to take the form 

d t F N = -M N F N + C®F N _ 2 (10) 

with the shorthand notation 

(C®F N - 2 )(t;r) = VC(^)F N _ 2 (f; ri , ,r N ) (11) 

n<m 

and with 

M * = -\Y, Va ^ r n-r m ) 8 rS d rt (12) 

n.m 

a/3 

being a symmetric (due to the incompressibility) , positive, singular elliptic 
differential operator of the 2 nd order. In the presence of the diffusion, Mjy 
should be modified by the subtraction of Xm K (^r n ) 2 - In the translation- 
ally invariant sector, i.e. for the homogeneous distributions of the scalar, 
the matrix T> a/3 (r) in the definition of Mn may be replaced by —D a @(r) 
since the contribution of V al3 (0) drops out. Below, we shall use this scaling 
form of dimension (length)^' 2 of the operators Mjy assuming also that the 
diffusion constant k has been taken to zero. 
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The Hopf equations (10) may be solved by induction with respect to N. 
Denoting by Pjv(i;r,r') the Green functions e~* Mjv (r, r'), we obtain 

F N (t;r) = J P N (t - t ; r,r') F N {t ; r') dr! 

+ [ ds [ P N (t-s; r,r') (C®F N - 2 )(s,r') dr! . (13) 
•'to 

For a concentrated initial distribution of 6 with fast decaying Appoint func- 
tions and for t —> oo, the Appoint functions F^(t;r) tend to the solution 
Fn (r) of the stationary version of the Hopf equations 

M N F N = C®F N _ 2 (14) 

inductively determined by the relations 

F N (r) = J G N (r,r') (C ® F N _ 2 )(r') dr! . (15) 

where Gjy(r, r') = J Q ^s Pn{ s ; £, r') is the kernel of the inverse of the oper- 
ator Mjy. In particular, the limiting stationary distribution is independent 
of the initial one and has vanishing odd-point functions. 



4. Zero mode dominance 

We are interested in the scaling properties of the stationary iV-point func- 
tions Fn(t) for the injection scale L large w.r.t. point separations, i.e. in 
the inertial range. Since Mjq have the scaling dimension of (length)^" 2 
and C{j} is approximately constant for large L, we might expect that the 
solutions of the chain (14) of equations scale like 

F N (Xr) = F N (r) . (16) 

This would be the normal Kolmogorov-Obukhov-Corrsin scaling. 

In the 1994 paper [3], Kraichnan has argued in favor of the anomalous 
scaling of the scalar correlators. His paper steered a renewed interest in 
the problem which led to the discovery of a simple mechanism by which 
the correlators avoid the normal scaling. As was realized in [4] and [5], see 
also [6] , for large L the normally scaling contributions to Fn are dominated 
by the ones of the zero modes of the operators Mjv- The stationary Hopf 
equations (14) leave the freedom to add the solutions of their homogeneous 
version M^fN = 0. It is not obvious, however, whether the solution (15), 
which is defined unambiguously, contains the scaling zero modes which 
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become leading in the inertial range. Such zero mode contributions were, 
indeed, found in the 2-point function. The latter may be calculated exactly 
[2] with the result 

F 2 (r) = A L 2 ^ + A ir 2 ~t + ... (17) 

where the terms vanishing as inverse powers when L — >• oo were omitted. 
The constant oc L 2 ~^ is a zero mode of operator M 2 and it dominates for 
r <^ L. The constant terms, as well as the ones independent of some of r n 's, 
however, do not contribute to the structure functions 



S N (r) = ((0(r)-0(O)n (18) 

of the scalar. As a result, the 2-point structure function scales as r 2 ~£, i.e. 
normally. 

What about the higher-point functions? In [4] and in [5] it was shown 
that the 4-point structure function is dominated for r <C L by the contri- 
bution of a zero mode of the operator M4 for, respectively, small velocity 
exponent £ and large space dimension d. In [7] and [8] these results have 
been extended to Appoint functions with the result 



S N (r) oc r Cjv with ( N = - < 



2{d+2) + )> ^ 



The scaling dimensions of the relevant zero modes where found by restrict- 
ing operators Mjy to scaling functions. Upon such a restriction, M^r's be- 
come operators with discrete spectrum to which standard perturbative tech- 
nics apply. A more difficult perturbative analysis around £ = 2 [6] [10] [11] 
confirms the zero mode dominance of the inertial-range scaling also in this 
regime. 



5. Slow modes 

The Green functions fW(i;r,r') = e _tMjv (r, r') have a natural interpreta- 
tion in the language of Lagrangian trajectories. They are the joint proba- 
bility distribution functions (PDF's) of the time t positions r n of the tra- 
jectories r n (s) starting at time zero at points r' n . This follows from the 
C = version of Eq. (13) if we recall that, in the absence of the source, 
the scalar density is carried along the Lagrangian trajectories. This in- 
terpretation of -P/v(t;r,r') shows that, effectively, the (differences of the) 
Lagrangian trajectories undergo a simple diffusion process with the space 
dependent diffusion matrix ^D a ^(r n — r m ) oc |r n — r m \^ so that they diffuse 
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very slowly when they are close but faster and faster when, eventually, they 
separate. 

Let us look closer at the effective diffusion of the Lagrangian trajectories. 
The scaling property of operators Mjy implies that 

P N (X 2 ^t; Ar,Ar ) = X~ Nd P N {t;r,r ) . (20) 

3±- 1 

Thus the distances scale as (time) as compared to (time) in the stan- 
dard diffusion. We may probe the effective spread of Lagrangian trajectories 
with scaling functions Jn such that /iv(Ar) = A CT /jv(r), e.g. with the sum 
of squares of distances between r n 's. The average of /at over the time zero 
positions of N Lagrangian trajectories which at time t pass by points r is 1 

(/jv) t , £ = J Pjv(t;r,r') Mr') dr' = \~° j P N (\ 2 ~H- Ar, r') f N (r') dr' 
= J Pn(t; C-)^r, r') f N (r') dr' (21) 

where we have used the scaling properties of Jn and Pn and have set 
r = \ 2 ~^t. If we let t grow and keep r constant, the last integral may be 
expected to tend to the limit / Pn(t; 0, r') /n(l') dr' so that, if the latter 
integral does not vanish, 

(f N } t , L oc C-)^ (22) 

This is the super-diffusive behavior which sets in for typical scaling probes 
/at. Suppose, however, that fx is a scaling zero mode of M^. Since 

J t (fN) t , r _ = ~ J PN(t;r,r!) (M N f N )(r!) dr' = 0, (23) 

the average (/7v)t,r is constant in time instead of growing like (-) 2 -«: the 
zero modes describe structures preserved in mean by the Lagrangian flow. 

A constant or linear functions are the obvious examples of scaling zero 
modes but a closer analysis [12] of operators shows that there is a 
whole discrete series /tv,i of them with scaling dimensions cij > 0. Besides, 
each such scaling zero mode gives rise to a tower of slow modes fN,ik, 
k = 0, 1, ... , starting at /at^o = fN,i, with scaling dimensions Oi + (2 — £)k. 
The averages (fN,ik)t,r are polynomials of order k in t so that they grow 

1 for the later convenience, we study the backward evolution in time 
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slower than the typical behavior (22). All these scaling modes appear in 
the small A asymptotic expansion [12] 

Pn(t; Ar,r') = ]T \°<+W k f N , ik (r) 5JV , ifc (r;r') (24) 

i,k 

which, when inserted on the right hand side of Eq. (21), gives 

(f"hr_ = E O^'" InMD J 9nMt;l') Ml') dr!. (25) 

i,k 

The functions gN,ik{ T ',L) are finite and decay fast for large r. For generic /at, 
the constant zero mode Jn,oo corresponding to gN,oo{T',L') = Pn{t; 0, r') 
gives the dominant contribution leading to the behavior (22). This term 
(and many others) vanishes for /jv = J'nji resulting in the slower growth. 

Upon integration over t, the asymptotic expansion (24) induces the one 
for the kernels of the inverses of Mtv's: 

Gn (Ar , r') = ]T f N>i (r) h Nii ( L ') (26) 

i 

(the contributions of the slow modes disappear under the time integration) . 
This expansion is behind the zero-mode dominance of the short-distance 
behavior of the iterative solutions (15) for the stationary AT-point functions 
of the scalar. The subtractions required in the passage from the correlation 
functions to the structure functions leave only the contributions of the zero 
modes /jv,i fully dependent on all r n 's. 

6. Fuzzy trajectories 

The asymptotic expansion (24) describes what happens when the final 
points of the Lagrangian trajectories are taken together. It was established 
with mathematical rigor for the case of two Lagrangian trajectories and was 
checked in perturbative approaches for AT > 2. At the first glance, it may 
look bizarre. Indeed, we may naively expect that if the final points r n of the 
trajectories converge then their initial points should also do with the joint 
PDF Pjv(t; Ar, r') becoming proportional to FJ <5( r n — r n+i) when A — > 0. 
Instead, Piv(t; Ar, r') tends to a regular limit gN,oo(t',L')- The mechanism 
by which P/v(f; r, r') avoids the singularity at r = is somewhat subtle [12]. 
Recall that the Lagrangian trajectories satisfy the differential equation (3). 
The uniqueness of solutions of such an equation requires the Lipschitz con- 
dition |v(i, r) — v(i, r')| ~ |r — r'| for r' tending to r. But our velocities are 
only Holder continuous in r with exponent | < 1, see Eq. (5). The result- 
ing non- uniqueness of the Lagrangian trajectories passing through a fixed 
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point violates the Newton-Leibniz paradigm and allows for a continuum 
of trajectories with coincident final points. Although the trajectories keep 
collapsing continuously, the probabilistic quantities such as the joint PDF's 
Pjv(t;r, r') still make perfect sense. They reflect, however, the fuzzyness of 
the trajectories in their asymptotics (24). 

A less rigorous-minded reader might object that the non-uniqueness of 
the Lagrangian trajectories is a mathematical nuisance since more realistic 
turbulent velocities, even when showing the behavior (5) with £ close to | in 
the inertial range, are smoothed on short distances by the viscous effects so 
that their Lagrangian trajectories are unique. A closer examination, shows, 
however, that such Lagrangian trajectories, although uniquely determined 
by sharp values at one time, exhibit a sensitive dependence on these val- 
ues signaled by non-zero Lyapunov exponents. As a result, if the fixed-time 
positions of the trajectories have an e-spread, the trajectories at different 
times are spread in a large region and will stay such if the viscosity is taken 
to zero first and the spread e only next. The non-uniqueness of the trajec- 
tories caused by their continuous collapse is simply a useful mathematical 
abstraction describing a real physical phenomenon: a fast spread of the tra- 
jectories in each realization of the high Reynolds number velocity field. Such 
a spread makes the identification of the trajectories practically impossible. 
It should be remarked that in the incompressible case, the collapse of the 
trajectories must be accompanied symmetrically by their branching since 
^V(^;r, r') = Piv(i;r' ,L)- I n compressible flows, this symmetry is broken. In 
particular, in the utmost compressible Burgers flows, the trajectories only 
collapse together (in a discrete process) sticking to the shocks. 

7. Numerical results 

The region of £ neither close to nor to 2 has up to now defied an exact 
analysis. The first numerical simulations of the system [13] [14] were based 
on the direct solution of the scalar equation (1). They did not cover the 
region of small £ and did not permit to discriminate between the small 
(cx £) anomalous structure-function exponents implied by the perturbative 
zero-mode analysis and the different (oc 1) predictions of a closure of the 
structure-function equations proposed in [3], see also [15]. 

Recently, a new numerical analysis [16] of the 4-point function of the 
scalar in two and three space dimensions has been performed for different 
values of £. It was based on a direct generation of samples of Lagrangian 
trajectories passing at time t by N fixed points. Such samples allow a 
direct simulation of the PDF's PN(t;r,r') as well as direct simulations of 
the stationary correlation functions F/v . Let us explain how the latter were 
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Figure 1. Anomalous exponent of the 4-point structure function of the scalar as a 
function of £ (U. Frisch, A. Mazzino and M. Vergassola, cond-mat/9802192) 



done in [16]. Upon setting to = — oo and 9(to) = in Eq. (4), one obtains 
for even N the relations 

N .t 

Fn(l) = (n/ ds f(s,r n (s))) = J2( \{T{v ni v m )) (27) 
«=i -°° n {n,ra}en 

where the sum, resulting from the average over the Gaussian source, is over 
the pairings II of the indices 1,...,N and where 

T(r n ,r m ) = f ds c C" (s) ' L r - is) ) (28) 

J — oo 

is approximately equal to the time spent within distance L by the two 

trajectories with the end-points r n and r m . The average of the product 

of such times may be simulated using the ensemble of trajectories passing 

through points r n . It is finite since the trajectories effectively separate as 
i 

(time) 2 -z . Due to the limited dependence of trajectories on the initial con- 
ditions discussed above, this remains true even if r„ and r m tend to each 
other. The comparison of -F/v(r) calculated this way for different configura- 
tions of points r n showed the dominance of the structure functions Sw(r) 
by subleading terms of Fjy(r), as in the discussion at the end of Section 5. 
The resulting values of the anomalous 4-point-structure-function exponent 
2C2 — C-4) presented in Fig. 1 borrowed from [16], confirm the predictions of 
the perturbative zero-mode analysis around £ = indicated in Fig. 1 by 
the broken-dotted line. 
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8. Conclusions 

We have exhibited a mechanism behind intermittency of the scalar struc- 
ture functions in the Kraichnan model of passive advection. The important 
point is that the anomalous scaling originated from a discrete series of solu- 
tions of the homogeneous Hopf equations without the forcing and diffusion 
terms. This, and the relation of such solutions to the effective fuzzyness of 
Lagrangian trajectories at high Reynolds numbers promise to be the stories 
that persist in other intermittent hydro dynamical systems. In particular, 
it has been indicated recently [17] how discrete solutions of the unforced 
Hopf equations may naturally give rise to a multiscaling picture of anoma- 
lous exponents in the Navier-Stokes case. 
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